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Southern rear edge populations of cold-adapted plants at low latitudinal margins will probably be critical 
for determining species’ responses to climate change. We focused on southern rear edge populations of 
an arctic-alpine plant, Dryas octopetala, in the high mountains of the Far East to assess genetic diversity 
and genetic structure by using allozyme variation and chloroplast DNA sequences. Allozyme data re¬ 
vealed low genetic diversity within populations (H c = 0.003-0.144), but high genetic differentiation 
among populations (G S t = 0.403), indicating characteristics of small and isolated populations that have 
persisted on mountain sky islands. The loss of genetic diversity was particularly pronounced in central 
Japan. On the basis of the cpDNA-haplotype distribution, the distinctive relationship between central 
Japan and not only northern Japan but also parts of Eurasia suggests that southernmost populations in the 
high mountains of central Japan are long-standing relicts of the Quaternary climate oscillations. Their 
refugial history therefore shaped the genetic diversity and genetic structure of the southern rear edge 
populations in Far East Asia. 
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Many arctic-alpine plants, so-called ‘circum¬ 
polar plants’, have vast ranges in the boreal belt of 
Eurasia and North America. Whereas the present 
distribution of arctic-alpine plants is beyond or 
above the potential tree line, fossil records indi¬ 
cate that arctic-alpine species spread southward 
into unglaciated areas to survive during Late- 
Quaternary glacial periods (Birks 2008). Since 
the last ice age, arctic-alpine plants have been mi¬ 
grating from their refugia upward to higher alti¬ 
tudes or northward to higher latitudes. Thus, the 


distribution of arctic-alpine plants has dramati¬ 
cally changed in response to Quaternary climate 
oscillations, causing them to experience frequent 
range expansions and fragmentations (Abbott & 
Brochmann 2003, Abbott 2008). In fact, many 
phylogeographic studies have revealed the genet¬ 
ic imprint of climate shifts in arctic-alpine spe¬ 
cies (e.g., Abbott et al. 2000, Skrede et al. 2006, 
Ikeda et al. 2008, 2009, Hirao et al. 2011, Winkler 
et al. 2012, Eidesen et al. 2013). 

The genetic imprint of range shifts are re- 
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markable in marginal populations. For arctic-al¬ 
pine species, populations on their northern mar¬ 
gin represent the expanding edges that colonize 
ice-free areas following glacial retreat (e.g., Alsos 
et al. 2007), whereas the southern marginal popu¬ 
lations represent the rear edges that have persist¬ 
ed or become extirpated (e.g., Chung et al. 2013). 
In contrast to the northern expanding popula¬ 
tions, which are newcomers, the southern rear 
populations are old and disproportionately signif¬ 
icant in terms of both evolution and conservation 
(Flampe & Petit 2005). Southern rear edge popu¬ 
lations of arctic-alpine species are isolated rem¬ 
nants of a more widespread distribution that ex¬ 
isted during the Quaternary cold stages, and are 
therefore called “glacial relicts” (Schmitt et al. 
2010, Flampe & Jump 2011). Southern relict pop¬ 
ulations have survived in the present interglacial 
refugia and even in some cases in long-standing 
refugia across glacial and interglacial periods. In 
addition, southern relict populations are typically 
small and isolated from the northern continuous, 
core populations. Their small size and prolonged 
isolation may possibly result in low levels of ge¬ 
netic diversity within a population and high lev¬ 
els of genetic differentiation among populations, 
thereby preserving high genetic distinctiveness. 
A survey of the population genetics of southern 
rear edge populations of arctic-alpine species 
may provide insights into the genetic conse¬ 
quences of historical isolation and fragmentation 
within plant species. 

In this study, intraspecific genetic diversity of 
southern rear edge populations of an arctic-alpine 
species, Dryas octopetala L. (Rosaceae), was in¬ 
vestigated to test the above hypotheses. Dryas 
octopetala is a long-lived shrub (life expectancy 
of more than 500 years; De Witte et al. 2012) and 
an important component of the arctic tundra and 
alpine heath in the Northern Hemisphere. Indi¬ 
viduals of D. octopetala form cushions by clonal 
proliferation. The flowers are self-fertile (Molau 
1993, Wada 1999), and the fruits are wind-dis¬ 
persed achenes with elongated, feathery styles. 
Dryas octopetala is widely distributed from the 
Eurasian and American arctic tundra to temper¬ 
ate mountains of Europe, North America, and 


Asia (Hulten 1959, 1968, Hulten & Fries 1986). 
Although Dryas octopetala shows greatest abun¬ 
dance at high latitudes, usually considered to be 
its core range, the southern populations at low lat¬ 
itudes are small and isolated. In a previous phylo- 
geographic study on D. octopetala sensu lato 
throughout most of its range, Skrede et al. (2006) 
found intrapopulation genetic diversity to be 
greater in northern populations than in southern 
populations. That study did not include samples 
from Far East Asia, however, where the intraspe¬ 
cific taxon, D. octopetala var. asiatica (Nakai) 
Nakai (Nakai 1932), occurs. The Asian endemic 
variety asiatica has more robust stems and short¬ 
er, broader, and blunter leaves than variety oc¬ 
topetala L. (Nakai 1932). Variety asiatica occurs 
in the high mountains of northeastern China 
(northern part of the Korean Peninsula) and in the 
Japanese Archipelago, where some of the south¬ 
ernmost populations are found. We therefore fo¬ 
cused on southern rear edge populations D. oc¬ 
topetala var. asiatica in the high mountains of 
Far East Asia to assess their genetic diversity and 
genetic structure. 

The unique phylogeographic history of arctic- 
alpine plants in Far East Asia compared with Eu¬ 
rope has been documented by, e.g., Ikeda et al. 
(2006). In the European Alps, arctic-alpine spe¬ 
cies survived glaciations during the Pleistocene 
ice ages either in nunatak refugia (e.g., Wester- 
gaard et al. 2011) or in peripheral-nonglaciated 
refugia (Brochmann et al. 2003, Alsos et al. 
2007). In contrast, Far East Asia remained free of 
a continuous cover of solid ice during the Pleisto¬ 
cene (Frenzel 1968, Svendsen et al. 2004), and 
high-altitude refugia have been crucial during 
warmer periods (interglacials) for survival 
against the upward encroachment of forest spe¬ 
cies. In fact, previous phylogeographic surveys 
revealed that the high mountains in the central re¬ 
gion of mainland Japan have been long-term re¬ 
fugia for a variety of alpine plant species (re¬ 
viewed in Fujii & Senni 2006) and have harbored 
unique haplotypes (e.g., Ikeda etal. 2008). There¬ 
fore, the unique refugial history should influence 
genetic diversity in rear edge populations of arc- 
tic-alpine species in Far East Asia. 
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Table 1. Overview of the collected materials of Dryas octopetala. 


ID 

Population 

Region 

Latitude/ 

Longitude 

Altitude 

(m) 

Year 

Collector(s)* 

Hi 

n 2 

Chroloplast DNA 
haplotype 

PI 

Warusawadake 

Central 

Honshu, Japan 

35°27’N/ 

138°09’E 

3000 

2005 

TM, MT 

30 

12 

A(12) 

P2 

Kisokomagatake 

Central 

Honshu, Japan 

35°47’N/ 

137°49’E 

2800 

2005 

NW, AH 

30 

10 

A(10) 

P3 

Yatsugatake 

Central 

Honshu, Japan 

35°58’N/ 

138°22’E 

2760 

2005 

NW, MM 

30 

10 

A(10) 

P4 

Suishodake 

Central 

Honshu, Japan 

36°25’N/ 

137°36’E 

2930 

2003 

NW, RM 

30 

- 

- 

P5 

Tateyama 

Central 

Honshu, Japan 

36°34’N/ 

137°37’E 

2700 

2003 

NW, RM 

60 

- 

- 

P6 

Yukikuradake 

Central 

Honshu, Japan 

36°47’N/ 

137°45’E 

2580 

2005 

NW, AH 

30 

10 

A(10) 

P7 

Taisetsuzan 

Hokkaido, 

Japan 

43°40’N/ 

142°55’E 

2060 

2006 

NW, KM 

30 

12 

B( 12) 

P8 

Changbaishan 

Northeast 

China 

42°02’N/ 

128°04’E 

2260 

2004 

NW, KK, QJL, 
XL 

30 

12 

C(12) 

P9 

Hemavan 

Scandinavian 
Peninsula, Sweden 

65°52’N/ 

15°14’E 

840-860 

2009 

AS 

- 

10 

F(3), G(5), H(l), 1(1) 


rii and n 2 are the number of individuals used in allozyme analysis and chroloplast DNA analysis, respectively. 

*Collectors: TM Takehiro Masuzawa, MT Miki Tomita, NW Naoya Wada, AH Akira Hirao, MM Madoka Miyazaki, RM Rie 
Matsuzato, KM Katsutoshi Matsunaga, KK Kunio Kawada, QJL Qi-Jing Liu, XL Xuefeng Li, AS Ayako Shimono. 



Fig. 1. (a) Sampling locations of the eight populations of Dry¬ 
as octopetala var. asiatica and a summary of cpDNA 
haplotypes. Population localities are summarized in Ta¬ 
ble 1. (b) A median-joining network of the seven cpDNA 
haplotypes detected in D. octopetala , including four 
haplotypes (F, G, H, I) from Sweden in northern Europe 
(see Table 1 for details). Small filled circles indicate 
missing haplotypes. Abbreviation: DD = Dryas drum- 
mondii. 

In this study, we surveyed the levels and dis¬ 
tribution of genetic diversity in Dryas octopetala 
var. asiatica, ranging in origin from northeastern 


China to the Japanese Archipelago, using 19 al¬ 
lozyme loci and sequences of two chloroplast 
DNA (cpDNA) regions. We tested our predictions 
of low genetic diversity and high genetic diver¬ 
gence in the southern rear edge populations on 
mid-latitude mountains of the Far East, consistent 
with the characteristics of climate relicts (Hampe 
& Jump 2011). 

Materials and Methods 

Sampling 

We sampled fresh leaves of Dryas octopetala var. 
asiatica from six populations in central Honshu 
(mainland Japan), including the southernmost 
population on Warusawadake, one population in 
Hokkaido (northern Japan), and one population in 
northeastern China, for both allozyme and 
cpDNA analyses (Table 1, Fig.la). All popula¬ 
tions were established on fellfields on wind¬ 
blown ridges where ericaceous dwarf shrubs, 
such as Vaccinium uliginosum L., were relatively 
common and where winter ground temperatures 
were approximately -15°C to -20°C (Wada et al. 
2006, Wada 2009). Fresh leaves from 30 to 60 
plants within each population were sampled. The 
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collections were from plants spaced at least 2 m 
apart to ensure selection from a broad area. The 
collected leaves were transported in a cooler and 
stored at -78°C prior to the experiments. 
Additionally, we obtained 10 specimens of Dryas 
octopetala var. octopetala from Hemavan, Swe¬ 
den, for cpDNA sequence analysis (see details in 
Table 1). The voucher specimens from Sweden 
were deposited in the herbarium of the Sugadaira 
Montane Research Center at the University of 
Tsukuba. Moreover, three samples of D. drum- 
mondii Rick, from Kananaskis, Canada, were 
used as an outgroup in the phylogenetic analyses. 

Allozyme variations 

The leaf tissue was homogenized in 1.0 mL of 
cold extraction buffer containing 0.1 mM Tris- 
HC1, 1 mM EDTA (4Na), 10 mM KC1, 10 mM 
MgCE, 0.4% 2-mercaptoethanol, and 10% poly- 
vinylpolypyrrolidone with the pH adjusted to 7.5. 
Enzymes were resolved on 7.5% polyacrylamide 
gels following the procedures of Shiraishi (1988). 
We analyzed the following 13 enzymes: aspartate 
aminotransferase (AAT; 2.6.1.1; three loci), acid 
phosphatase (ACP; 3.1.3.2; one locus), alcohol de¬ 
hydrogenase (ADH; 1.1.1.1.; two loci), diaphorase 
(D1A; 1.6.4.1; two loci), glucose-6-phosphate de¬ 
hydrogenase (G6PD; 1.1.1.49; one locus), isoci¬ 
trate dehydrogenase (IDH; 1.1.1.42; one locus), 
leucine aminopeptidase (LAP; 3.4.11.1; two loci), 
malate dehydrogenase (MDH; 1.1.1.37; one lo¬ 
cus), phosphoglucose isomerase (PGI; 5.3.1.9; 
one locus), phosphoglucomutase (PGM; 5.4.2.2; 
two loci), 6-phosphogluconate dehydrogenase 
(6PGD; 1.1.1.44; one locus), shikimate dehydro¬ 
genase (SKD; 1.1.1.25; one locus), and triose- 
phosphate isomerase (TPI; 5.3.1.1; one locus). We 
resolved the putative loci under the assumption of 
Mendelian inheritance because Dryas octopetala 
sensu lato is uniformly diploid with 2 n = 18 (Elk- 
ington 1971). 

The allele frequencies in each population 
were calculated from the 19 loci encoding the 13 
enzyme systems. A locus was considered poly¬ 
morphic if the most common allele had a fre¬ 
quency less than 0.99. We calculated standard 
measures of genetic variation, including the per¬ 


centage of polymorphic loci (P), mean number of 
alleles per locus (A), and expected heterozygosity 
or gene diversity (II C ). Deviations from Hardy- 
Weinberg expectations were quantified for each 
polymorphic locus within each population with 
Wright’s fixation index ( F ). Fixation indices were 
tested for significance with the chi-squared anal¬ 
ysis (Li & Horvitz 1953). The amount of genetic 
differentiation among populations at the gene 
pool level was estimated using the G SJ measure 
(Nei 1973) in which gene diversity across all pop¬ 
ulations ( H t ) was calculated for each locus and 
then averaged over all loci (Culley et al. 2002). To 
visualize similarities between populations, we 
performed a principal coordinate analysis (PCO) 
based on Nei’s standard genetic distance with 
sample size bias correction (Nei 1978) using the 
R software package (R Core Team 2013). The re¬ 
lationships between genetic and geographic dis¬ 
tance (isolation by distance) among all popula¬ 
tions were estimated (Slatkin 1987, 1993, Rousset 
1997) by correlating G S t/(1-G S t) with the natural 
logarithm of geographic distance (km) in a Man¬ 
tel test with 9999 permutations using the R soft¬ 
ware package. 

Sequence of chloroplast DNA 

Two non-coding cpDNA regions were se¬ 
quenced: the atpB-rbcL spacer (Hodges & Ar¬ 
nold 1994) and the trnV-ndhC spacer (Shaw et al. 
2007). We analyzed 10 to 12 plants per popula¬ 
tion, with the exception of the Tateyama and Su- 
ishodake populations, from which every sample 
was used for the prior allozyme analysis. DNA 
was extracted by the cetyltrimethyl-ammonium 
bromide (CTAB) method (Stewart Jr & Via 1993) 
or by using the DNeasy Plant Mini Kit (Qiagen, 
Hilden, Germany) following the manufacturer’s 
protocol. PCR amplification was performed in a 
10-pL reaction mixture containing 0.5 pL of tem¬ 
plate DNA, 5 pi of 2* AmpliTaq Gold 360 Master 
Mix (Applied Biosystems, Foster City, CA, USA), 
and each pair of primers at 0.2 pM: atpB and rbcL 
(Terachi 1993), and tr«F (UAU x2 and ndhC (Shaw 
et al. 2007). The amplification parameters were 
95°C for 10 min; 35 cycles at 95°C for 1 min, 
50°C for 1 min with a ramp of 0.3°C/s, and 65°C 
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for 4 min; and 65°C for 5 min following the 
“ rpll6 ” cycling program (Shaw etal. 2005). The 
PCR products were purified using an ExoSAP-lT 
Kit O.lx (GE Healthcare, Little Chalfont, UK). 
The cycle sequencing reactions of forward and 
reverse strands were performed with the same 
PCR primers, using O.lx diluted BigDye v.3.1 
(Applied Biosystems). The sequencing products 
were analyzed on an AB1 Prism 3130 automated 
sequencer (Applied Biosystems). Forward and re¬ 
verse sequences were assembled using BioLign 
(Hall 2001). The DNA sequences were aligned 
using MEGA6 (Tamura et al. 2013) and the mul¬ 
tiple sequence alignment was refined manually. 
Indels were treated as single characters resulting 
from one mutation event. Haplotype was deter¬ 
mined based on a combination of the two non¬ 
coding cpDNA regions. Haplotype diversity (/?) 
and nucleotide diversity (n) were calculated at the 
population level using the program DnaSP ver¬ 
sion 5.0 (Librado & Rozas 2009). A median-join¬ 
ing network with Dryas drummondii as an out¬ 
group was inferred using NETWORK v. 4.6.0.0 
(Bandelt et al. 1999), with the indels coded as a 
fifth state. 

Results 

Allozyme diversity 

Nineteen loci (AAT-1, AAT-2, AAT-3, ACP, 
ADH-1, ADH-2, DIA-1, DIA-2, G6PD, IDH, LAP- 
1, LAP-2, MDH, PG1, PGM-1, PGM-2, 6PGD, 
SKD, and TPI), coding for 13 enzymes, were re¬ 
solved (Appendix I). Seven loci (ACP, DIA-2, 
IDH, LAP-1, MDH, PGM-1 , and TPI) were mono- 
morphic, whereas 12 loci ( AAT-1, AAT-2, AAT-3, 
ADH-1, ADH-2, DIA-1, G6PD, LAP-2, PGI, 
PGM-2, 6PGD, and SKD) showed considerable 
polymorphism in at least one of the populations. 
Most alleles were shared among populations, but 
12 private alleles (sensu Slatkin 1985) were pres¬ 
ent (Table 2, Appendix I). Ten of the 12 private 
alleles were detected on Changbaishan, one was 
detected on Taisetsuzan, and another one was de¬ 
tected on Yatsugatake. 

Southern rear edge populations of Dryas oc¬ 
topetala var. asiatica showed low levels of allo¬ 


zyme diversity that varied substantially among 
individual populations (Table 2). The mean per¬ 
centage of polymorphic loci (P) across all popu¬ 
lations was 0.684, with values at the population 
level ranging from 0.105 to 0.579 with a mean of 
0.224. Overall, the mean number of alleles per lo¬ 
cus (A) was 1.9, and the estimates at the popula¬ 
tion level ranged from 1.1 to 1.7 with a mean of 
1.2. The observed heterozygosity ( H 0 ) ranged 
from 0.011 to 0.160 with a mean of 0.042. Al¬ 
though locus-specific significant deviations from 
Hardy-Weinberg expectations were found (four 
positive, two negative) among 26 tests, the mean 
fixation index ( F ) of each population was not sig¬ 
nificant (Appendix II). The gene diversity across 
all populations (H T ) was 0.070, whereas the val¬ 
ues at the population level (H e ) ranged from 0.003 
to 0.144 with a mean of 0.042. The population at 
Changbaishan in northeastern China had the 
highest level of genetic diversity (0.144) while 
populations in the Japanese Archipelago exhibit¬ 
ed H e values (0.003-0.068) less than half of that 
(Table 2). Within the populations from the Japa¬ 
nese Archipelago, the population on Yatsugatake 
(P3) had the highest level of genetic diversity. 
The estimated genetic differentiation among pop¬ 
ulations, G st , was 0.403, indicating that approxi¬ 
mately 40% of the total genetic diversity resulted 
from differences among populations. 

The genetic similarities between populations 
were analyzed by means of PCO (Fig. 2). The first 
component (PCOl), which accounted for 73.2% 
of the variation, separated the populations in Ja¬ 
pan (P1-P7) from the Changbaishan population 
in northeast China (P8). The second component 
(PC02), which accounted for 14.1% of the vari¬ 
ance, separated the Taisetsuzan population in 
northern Japan (P7) from the remainder of the 
populations. In the scatter plot, we identified a 
cluster of six populations in central Japan (Pl- 
P6) that were situated between the northern Ja¬ 
pan population (P7) and the northeastern China 
population (P8). Within populations in central Ja¬ 
pan (P1-P6), however, the Yatsugatake popula¬ 
tion (P3) was plotted at a distance from the re¬ 
maining five populations. Population differentia¬ 
tion followed a general pattern of isolation by dis- 
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Table 2. Allozyme diversity of the eight populations of Dryas octopetala var. asiatica in Far East Asia. 


ID 

Population 

P 

A 

A 

-^private 

Ho 

H c 

PI 

Warusawadake 

0.211 

1.2 

0 

0.047 

0.048 

P2 

Kisokomagatake 

0.158 

1.1 

0 

0.021 

0.026 

P3 

Yatsugatake 

0.316 

1.3 

1 

0.053 

0.068 

P4 

Suishodake 

0.105 

1.1 

0 

0.004 

0.003 

P5 

Tateyama 

0.158 

1.1 

0 

0.013 

0.012 

P6 

Yukikuradake 

0.105 

1.1 

0 

0.011 

0.009 

P7 

Taisetsuzan 

0.158 

1.1 

1 

0.025 

0.024 

P8 

Changbaishan 

0.579 

1.7 

10 

0.160 

0.144 


Average 

0.224 

1.2 

1.5 

0.042 

0.042 


Overall 

0.684 

1.9 

12 

0.042 

0.070 


P is proportion of polymorphic loci, A is the mean number of allele per locus, A piivate is the number of “private allele” (Slat- 
kin!985), H 0 is the observed heterozygosity, and He is the expected heterozygosity or gene diversity. 



PC01(73.2%) 


Fig. 2. Principal coordinate analysis of the eight populations 
of Dryas octopetala var. asiatica based on Nei’s genetic 
distance obtained from allozyme analysis. Circles, the 
triangle, and the square show the populations in central 
Japan (P1-P6), the population in northern Japan (P7), 
and the population in northeast China (P8), respectively. 
See Table 1 and Figure 1 for details regarding the popu¬ 
lations. 


tance; Mantel test showed a significant pattern 
over the entire geographic range (r = 0.597, p < 
0.01: Fig. 3). 

cpDNA haplotypes and their geographic distri¬ 
bution 

The concatenated sequence of the two non¬ 
coding regions of cpDNA had a length of 1173 
bp, of which 13 sites were variable with five sub¬ 



Geographic distance (km) 

Fig. 3. Pattern of isolation by distance for all pairs of the eight 
populations of Dryas octopetala var. asiatica. 


stitutions and eight indels. Three haplotypes (A, 
B, and C) were identified across 64 samples of 
Dryas octopetala var. asiatica from Far East 
Asia (Table 1, Fig. 1). Additionally, four hap¬ 
lotypes (F, G, H, and I) were detected in sam¬ 
ples from Sweden in northern Europe. The se¬ 
quences are available at DDBJ/EMBL/GenBank 
(AB685700-AB685705, AB720971, AB720972, 
and AB720975-AB720977). The median-join¬ 
ing network revealed three distinct clades (Fig. 
lb). One specific clade contained only the single 
haplotype (A) from central Japan. The second 
clade, comprising haplotypes B and C, occurred 
in northern Japan and northeast China, respec¬ 
tively. Finally, an additional clade, which includ- 
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ed haplotype F, G, H, and I, was found in Sweden. 
Haplotype F was located in the central part of the 
network. All populations in Far East Asia were 
fixed by a single haplotype, whereas the popula¬ 
tion in Sweden was genetically divergent (haplo¬ 
type diversity, h = 0.711; nucleotide diversity, n = 
0.00048). 

Discussion 

The geographic distribution of the cpDNA 
haplotypes in Far East Asia indicates an exten¬ 
sive genetic structure among central Japan, 
northern Japan, and northeast China, in each of 
which a unique haplotype was fixed. All popula¬ 
tions in central Japan were of cpDNA haplotype 
A, which is phylogenetically distinct not only 
from haplotype B in northern Japan but also from 
haplotype C in northeast China (Fig. 1). This re¬ 
sult is consistent with previous phylogeographic 
studies of Japanese alpine plants, which hypoth¬ 
esized that the high mountains of central Japan 
served as long-standing isolated refugia during 
Quaternary climatic oscillations (reviewed in Fu- 
jii & Senni 2006). Thus, the southernmost popu¬ 
lations of central Japan should be older relicts 
that have survived more than one glacial/ inter¬ 
glacial cycle. Altitudinal shifts in mountains over 
3000 m a.s.l. in central Japan may allow long¬ 
term population persistence of arctic-alpine 
plants, some of which have speciated in situ (Ike- 
da et al. 2012, 2014). In Far East Asia, moreover, 
the lack of extensive glaciation (Frenzel 1968, 
Svendsen et al. 2004) has contributed to the pres¬ 
ence of long-standing refugia. In Europe, howev¬ 
er, southern alpine populations have a different 
evolutionary history because glaciers obliterated 
most of their previous areas of distribution. Al¬ 
though some phylogeographic studies have re¬ 
ferred to glacial relicts of arctic-alpine species in 
the Alps (e.g., Reisch et al. 2003), a survey of the 
literature revealed that most phylogeographic 
patterns of arctic-alpine species in Europe are the 
result of recent (in most cases postglacial) vicari- 
ance or dispersal events (Schmitt et al. 2010). 
Therefore, the long-standing refugia in central 
Japan provide a unique phylogeographic history 


along with the impetus for shaping the genetic di¬ 
versity and genetic structure of the southern rear 
edge populations in Far East Asia. 

Allozyme data revealed low genetic diversity 
in Dryas octopetala var. asiatica within the over¬ 
all population ( H T = 0.070) and at the population 
level (H c = 0.003-0.144) (Table 2). Of particular 
note is that the level of genetic diversity in south¬ 
ern populations was much lower than has been 
reported for more northern populations around 
the North Atlantic (H c = 0.057-0.423; Philipp & 
Siegismund 2003) and for other boreal-temperate 
and woody plants (H e = 0.272 and 0.298; Hamrick 
& Godt 1989). The low genetic diversity in D. oc¬ 
topetala var. asiatica is supported by the cpDNA 
data, although based on limited sampling (h = 0, 
n= 0 in Asian Far East populations, and h = 0.711, 
7r = 0.00048 in a Sweden population). 

In contrast to the low level of genetic diver¬ 
sity within populations, high genetic differentia¬ 
tion was found among populations of Far East 
Asia (G st = 0.403). As supported by the cpDNA 
findings, allozyme divergence among the popula¬ 
tions has increased in rear edge populations in the 
Far East, in contrast to northern populations of 
the same species: F ST = 0.00-0.11 in Alaska (Max 
et al. 1999); r/) ST (F ST analogue by AMOVA) = 
0.216 in countries around the North Atlantic (Vik 
et al. 2010). Additionally, we found some evi¬ 
dence of isolation by distance based on allozyme 
analysis (Fig. 2), indicating that historical gene 
flow has been restricted by geographical distance 
but was not obliterated across Far East Asia, 
where topographical barriers provided by moun¬ 
tains and/or the sea may explain the limitation of 
gene flow. The southern rear edge populations of 
Dryas octopetala are located on montane sky is¬ 
lands on true islands in the Far East. The spatial 
isolation of populations inhabiting mountain sky 
islands enforces genetic differentiation among 
populations (DeChaine & Martin 2005). Geo¬ 
graphically, northeastern China is part of the 
Eurasian continent, whereas the Japanese Archi¬ 
pelago is an island chain located along the East 
Asian littoral region, where populations are hy¬ 
pothesized to have been isolated from those on 
the continent. The PCO analysis, in which the 
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first component separated the Japanese popula¬ 
tions from the Chinese population, and the sec¬ 
ond component, which separated northern Japan 
from central Japan, supports isolation by sea. 
Geographic features serving as barriers thereby 
accelerate genetic differentiation among popula¬ 
tions. 

The low levels of genetic diversity were par¬ 
ticularly pronounced in the populations of central 
Japan, whereas the population on Changbaishan 
in northeast China retained considerable genetic 
diversity (Table 2). The levels of within-popula- 
tion genetic diversity appear to be correlated with 
population size. Most of the populations in cen¬ 
tral Japan are small. For example, approximately 
150 patches of Dryas octopetala var. asiatica oc¬ 
cur on Tateyama (Fig. 1) (Wada 2008). The popu¬ 
lations on Yatsugatake are exceptional in com¬ 
prising relatively many patches along the main 
ridge from Akadake to Yokodake (A. S. Hirao & 
N. Wada, personal observation) and have a higher 
level of genetic diversity than do other popula¬ 
tions in central Japan (Table 2). On Changbaishan 
in northeast China, where Dryas octopetala var. 
asiatica occurs from the sub-alpine to the alpine 
zone, the population is considerably larger. Rela¬ 
tively highly polymorphic loci and most of the 
private alleles were found in population from 
northeast China. Additional samples, especially 
outside of the Japanese Archipelago, are neces¬ 
sary to further clarify the trends in genetic diver¬ 
sity of the Asian endemic variety asiatica. 

The reproductive system of plants has long 
been hypothesized to influence their genetic di¬ 
versity (Hamrick & Godt 1989, Nybom & Bartish 
2000). In fact, the difference in genetic diversity 
in two arctic-alpine plants, Diapensia lapponica 
L. and Empetrum nigrum L., at the southernmost 
edge of their range is mainly attributable to dif¬ 
ferent life histories and reproductive traits (Chung 
et al. 2013). Thus, the reproductive system of D. 
octopetala, such as self-fertility (Molau 1993, 
Wada 1999), may also affect genetic variation. 
Accelerated inbreeding with decreasing popula¬ 
tion size leads to an increase in homozygocity, 
genetic drift, and random changes in allele fre¬ 
quencies. In turn, if the Japanese populations self 


pollinate and the Changbaishan populations out- 
cross, trends in the level of genetic diversity may 
occur. We did not observe the predicted pattern in 
the inbreeding coefficient (Appendix II), howev¬ 
er, possibly because of the extremely low poly¬ 
morphism within the Japanese populations in 
which most alleles are fixed. Application of high¬ 
ly polymorphic markers, such as microsatellites, 
can identify the relationship between inbreeding 
and genetic diversity in southern rear edge popu¬ 
lations. 

Changes in climate threaten many species, 
particularly those in arctic-alpine regions (IPCC 
2007, Engler et al. 2011, Pauli et al. 2012). Even 
within arctic-alpine species with wide-ranging 
distributions, there is a high risk of losing south¬ 
ern populations due to climate warming (Thuiller 
et al. 2005). In fact, the fossil record shows that 
rear edge populations are often extirpated (Davis 
& Shaw 2001). Furthermore, the southern popu¬ 
lations may be genetically vulnerable, thus their 
viability may be overestimated. Low levels of ge¬ 
netic diversity at the rear edge of species’ ranges 
may limit evolutionary potential in the context of 
climate change. Further clarification of the evolu¬ 
tionary history of Dryas octopetala var. asiatica 
can improve our understanding of the complicat¬ 
ed response of arctic-alpine species to historical 
and contemporary climate change. 
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